Climate-smart rice (Oryza sativa L.) genotypes identification using stability analysis, multi-trait selection index, and genotype-environment interaction at different irrigation regimes with adaptation to universal warming

Climate change has brought an alarming situation in the scarcity of fresh water for irrigation due to the present global water crisis, climate variability, drought, increasing demands of water from the industrial sectors, and contamination of water resources. Accurately evaluating the potential of future rice genotypes in large-scale, multi-environment experiments may be challenging. A key component of the accurate assessment is the examination of stability in growth contexts and genotype-environment interaction. Using a split-plot design with three replications, the study was carried out in nine locations with five genotypes under continuous flooding (CF) and alternate wet and dry (AWD) conditions. Utilizing the web-based warehouse inventory search tool (WIST), the water status was determined. To evaluate yield performance for stability and adaptability, AMMI and GGE biplots were used. The genotypes clearly reacted inversely to the various environments, and substantial interactions were identified. Out of all the environments, G3 (BRRI dhan29) had the greatest grain production, whereas G2 (Binadhan-8) had the lowest. The range between the greatest and lowest mean values of rice grain output (4.95 to 4.62 t ha-1) was consistent across five distinct rice genotypes. The genotype means varied from 5.03 to 4.73 t ha-1 depending on the environment. In AWD, all genotypes out performed in the CF system. With just a little interaction effect, the score was almost zero for several genotypes (E1, E2, E6, and E7 for the AWD technique, and E5, E6, E8, and E9 for the CF method) because they performed better in particular settings. The GGE biplot provided more evidence in support of the AMMI study results. The study's findings made it clear that the AMMI model provides a substantial amount of information when evaluating varietal performance across many environments. Out of the five accessions that were analyzed, one was found to be top-ranking by the multi-trait genotype ideotype distance index, meaning that it may be investigated for validation stability measures. The study's findings provide helpful information on the variety selection for the settings in which BRRI dhan47 and BRRI dhan29, respectively, performed effectively in AWD and CF systems. Plant breeders might use this knowledge to choose newer kinds and to design breeding initiatives. In conclusion, intermittent irrigation could be an effective adaptation technique for simultaneously saving water and mitigating GHG while maintaining high rice grain yields in rice cultivation systems.


The variance in rice grain yield
This finding is supported by the range in the mean yield of the genotypes that were evaluated, which varied for AWD from 3.8 (t ha -1 ) (G2 in environment E3) to 5.4 (t ha -1 ) (G3 in environments E2 and E7, G4 & G5 in environments E1) (Fig. 1A,C).The G2 genotype was discovered to have the lowest yield at one place, E3, whereas the G4 genotype was detected at locations in E7, E1, E2, and E4 (Fig. 1).In every area, genotype G3 produced the most (Fig. 1A).The observed variation in the mean yield of the examined genotypes for CF (Fig. 1B,D) ranges from 3.9 (t ha -1 ) (G5 in environment E6) to 5.4 (t ha -1 ) (G3 in environments E2, E8, and E9).In most of the sites, the G3 genotype yield in both water regime treatments is consistent (Fig. 1A-D).

Values of the combined and AMMI model's analysis of variance (ANOVA) applied to yield and yield-related characteristics
The combined and AMMI model's Analysis of Variance (ANOVA) for AWD and CF systems are shown in Tables 1 and 2. The results of the combined analysis of variance showed that the genotypes, environments and the GEI were significant (Tables 1 and 2).The current study used the analysis of the variance function of the Figure 1.The rice grain yield (t ha -1 ) performance variation of studied five (5) rice genotypes across nine (9)  environments in both AWD = Alternate wet and dry (A, C) and CF = continuous flooding (B, D) conditions.AMMI model to analyze the yield and yield performance (PH, ET, PL, TSW, DM, and GY) of five prospective rice genotypes evaluated in nine environments throughout the dry season.The ANOVA showed significant impacts from the environment, G × E interactions, and genotype for each attribute.Based on the AMMI analysis of variance, every characteristic in systems PC1 and PC2 was determined to be significant.There was variation in the genotype rankings between settings for each feature that was being studied.Given that GEI has a considerable impact on each of the studied characteristics except for ET and PL; AMMI's model was employed in conjunction with multiplicative effects analysis to identify stable genotypes.The AMMI model with two significant IPCs is the best-predicted model 43,44 .At the 1% probability level, the first three components of PH were significant and explained 99.5% and 93.5% of the changes related to GEI for the AWD and CF water regimes, respectively.The third factor was the effective tiller (ET), which accounted for 95.8% and 94.8% of the GEI impact in the AWD and CF scenarios, respectively.Significant findings were found for each explanation at the 1% probability level.A significant component, or 89.7% and 93.9% for AWD and CF, respectively, supported the PL-related variations at the 1% probability level.The two main components of the TSW explained 89.8% and 81.2% of the GEI impact at the AWD and CF levels, respectively.In terms of DM, the GEI for the AWD and CF water regimes is Table 1.Combined ANOVA and AMMI ANOVA applied to yield and yield-related trait values of rice genotypes tested at AWD conditions within nine (9) environments.GEN = Genotype; ENV = Environment; MS = Mean Square; AMMI = Additive main effects and multiplicative interaction; GEI = Genotype by environment interaction; PC = Principal component; *** = Significant at 0.1%;** = Significant at 1%; ; ** = Significant at 1%; * = Significant at 5% probability level; df = Degrees of freedom.www.nature.com/scientificreports/accounted for by the first three interaction principal components of AMMI and three PCs with highly significant differences (p < 0.001), respectively.In this study, PC1 (59.8% and 66.6%) and PC2 (25.7% and 25%) each made a substantial contribution to the understanding of the G × E interaction for GY.The sums of squares for PC1 and PC2 together account for 85.5% and 91.6% of the total GEI at both the AWD and CF levels.

Genotypes and water regimes are graphically shown in the AMMI Biplots
According to the research, the distribution of genotypes and environments is represented by the AMMI1 biplot stability of yield, where PC1 accounts for 66.6% in CF systems and 59.1% in the AWD approach (Fig. 2A,B).To determine the mega-environment and stable genotypes throughout AWD and CF, biplots for AMMI1-Means vs. PC1 and AMMI2-PC1 vs. PC2 were created (Figs.2A,B).The means of the environments and genotypes were plotted against their IPCA1 using the AMMI1 biplot.The low-yielding nature of environments E6, E3, and E5 was later shown (Fig. 2A AMMI1 Biplot).The genotypic and environmental scores that correspond to the GY IPCA1 and IPCA2 are shown in the AMMI2 biplot (Fig. 2C,D AMMI2 Biplot).The AMMI 1 biplot is typically interpreted as follows: displacements along the ordinate indicate differences in interaction effects, while displacements along the abscissa reflect variations in main (additive) effects.Zero PC1 genotypes are more adaptable to all environments and are less influenced by external variables.As evidenced by their nearly low PC1 scores, varieties G5 were the most stable and adaptive genotypes under these circumstances (Fig. 2A).However, since genotype G3's mean yield was higher than the other genotypes' , it is recommended.In conclusion, the best yield might not come from a sturdy and flexible variety.The first two principal components, PC1 and PC2, could account for 84.8% of the total variance in the GEI (Fig. 2C).Regarding the grain yield attribute, there was the least amount of interaction between genotypes and environments in G1 > and G2.The AMMI1 biplot (Fig. 2B) showed that grain yields were below average for G3, G2, and G5, but they were above average for G1 and G5.G5 was the genotype that produced the least from these.Strong G × E interaction effects were seen in G5 and G3 across the genotypes under evaluation.Conversely, G1 displayed a reduced G × E interaction.The grain yield display interaction of PC1 and PC2 of the genotypes that were assessed in the eight settings is represented by the AMMI2 biplot in Fig. 2D.Yan and Kang 45 suggest that distances from the biplot origin are a useful measure of the degree of genotype or environment interaction.While they did not exert considerable interaction forces, the environmental vectors with short spokes in E1, E9, E3, E4, and E5 in this figure greatly contributed to the genotype's stability.The environmental vectors with long spokes, on the other hand, show high interaction.As a result, environments E7 and E8 had long spokes, demonstrating their high discriminating ability.

The examination of rice grain yield stability
Every genotype in the CF and AWD systems performed very identically in each of the nine settings (Tables 3  and 4).Table 3 shows that, in all circumstances, G3 and G2 had the highest and lowest grain yields, respectively.The greatest and lowest mean values (4.95 to 4.62 t ha -1 ) of the variation in rice grain yields across the five rice genotypes were similar to one another.Table 4 displays the genotype means, which varied between 5.03 and 4.73 t ha -1 depending on the environment.Tables 3 and 4 show that all genotypes performed better in AWD than in the CF system.Using the AMMI biplot adaptation map (Nominal plot) for genotypes or generic genotypic adaptation.According to the AMMI biplot adaption map of yield, the AWD approach was more suited for G1 (Binadhan-10) and G3 (BRRI dhan-29).These two genotypes exhibited comparable performance in all situations.Although it varies from environment to environment, G5 (Binadhan-17) performed best in environment-3 (Union: Kolma).Only the G1 (Binadhan-10) exhibited reduced G × E interaction in the CF system.Similar to the AWD approach, Binadhan-17 also had the largest G × E interaction in this instance (Fig. 3).

Using the which-won-where method of the polygon of GGE biplot methodology, the specific genotypic adaption
The vertex genotype in the AWD approach was G4 (BRRI dhan47), while the genotype with the highest responsiveness was G2 (Binadhan-8), per the study's Polygon of the GGE biplot of rice grain yield.In the CF system, G5 (Binadhan-17) has been found in both responsive and vertex genotypes.The genotype main effects and the GE interaction, which were 54.94% and 56.31%, respectively, were explained by the PC1 in the AWD technique and in the CF system (Fig. 4).

GGE biplot genotype view with the average environment coordination (AEC) of rice grain yield
Figure 5A,B show the genotype ranking according to mean performance and stability.According to previous research, PC2 (a measure of instability) of a GGE biplot must approximate the GE effects associated with each genotype if PC1 of the biplot approximates the genotype's main effects 34,46 .The average environment coordinate (AEC) axis, which is defined by the average PC1 and PC2 scores of all the environments, is the line that goes through the biplot origin and the average environment.The yield stability of genotypes was assessed using the AEC approach, which uses the average principal components in all conditions.Genotype G3 (BRRI dhan29) exhibited the highest grain yield and excellent stability under Alternate Wetting and Drying (AWD) conditions, as depicted in Fig. 5A.Conversely, genotype G1 (Binadhan-10) demonstrated average grain yield but maximum stability.Genotypes G4 (BRRI dhan47) and G5 (Binadhan-17) showed grain yields close to that of G3 (BRRI dhan29) but with lower stability.Overall, it is inferred that genotypes G1 (Binadhan-10) and G3 (BRRI dhan29), offering comparable grain yields to optimal genotypes while maintaining sufficient stability, are the most desirable.Concerning the continuous flooding system (CF), genotype G3 (BRRI dhan29) displayed a nearly average grain yield, while G4 (BRRI dhan47) exhibited a high grain yield with lower stability.In contrast, G1 (Binadhan-10) showed average grain yield with the highest stability.Genotypes G2 (Binadhan-8) and G5 (Binadhan-17) were associated with the lowest grain yield and stability.In conclusion, Binadhan-10 and BRRI dhan29 emerge as the superior genotypes, demonstrating grain yield comparable to the optimal genotype while exhibiting adequate stability across conditions.

Analyzing genotypes in relation to the ideal genotype
The GGE biplot model is an interesting tool for evaluating genotypes against an ideal genotype.Concentric circles were added to a genotype-focused scaling GGE biplot graph 35,47 to more clearly illustrate the difference between genotypes and the ideal genotype.The GGE biplot graph (Fig. 6A,B) demonstrates that genotype G3 was positioned first in the concentric circle in both the AWD and CF techniques.Therefore, G3 was the best genotype position-better than any other rice genotype-and was followed by G4, G1, and G3.The genotypeenvironment interaction shown in Fig. 6A and the GGE biplot analysis of rice grain yield indicate that G3 (BRRI dhan-29), which is nearly the ideal genotype theoretically, is the best genotype.They were therefore regarded as being outside of Fig. 6B's ideal inbred group.G5 (Binadhan-17) in the CF system and the inbred G2 (Binadhan-8) in the AWD technique, on the other hand, were found to be extremely distant from the optimal genotype.

Assessment of environment in comparison to ideal environment
In optimal conditions, the longest vector-which had a low IPCA-was discovered to lie inside the centers of concentric circles.The environment-focused GGE biplot's second concentric circle depicts the ideal environment, and the settings that are most like it are the ones that are desired.In this sense, E2 is in the second concentric circle and has been in the ideal environment (Fig. 7A).Though they may be viewed as unfavorable (having fewer examples) for selecting cultivars that are widely adapted, the E6 and E9 environments were far from ideal conditions, thus selecting particularly suitable cultivars may benefit from them.Despite not being statistically representative, the most discriminating samples in the current analysis were E1, E2, E6, and E7 from Fig. 7A.Their long projection against the ATC Y-axis and distance from the plot source for the AWD approach were probably the causes of this.Consequently, E9 and E5, which represent the ideal conditions based on this research, are shown in Fig. 7B.An environment that is most representative of all environments and most effective at differentiating between genotypes is considered optimal.Similarly, because they were closer to the ideal environment, E4, E2, and E8 were judged to be the second strongest in terms of genotype discrimination (Fig. 7B).On the other hand, environment E7 was found to be distant from the ideal environment and was believed to be less successful in distinguishing between genotypes (Fig. 7B).On the E9, the standard environment is displayed.Figure 7B demonstrated that while E5, E6, E8, and E9 were not statistically representative in the current investigation, they were the most discriminating because they showed the longest projection over the ATC Y-axis and were located farthest from the plot source for the CF technique.

Correlation between various stability statistics
To determine the relationship between each pair of AMMI stability parameters, Spearman's correlations (Fig. 8) discovered a strong association between ASV and sij, Ecoval and sij, ASV and Ecoval, HMRPGV, and RPGV.In both AWD and CF systems, ASV, sij, HMRPGV, and RPGV exhibited a strong correlation with each other, while the other measures showed a non-significant correlation.The results showed that in both the CF and AWD systems, the same parameters showed a significant correlation in the same patterns.A heatmap-based Speraman's correlation coefficient was used to analyze the relationships between mean grain yield and the calculated stability indicators (Fig. 8).The results demonstrated a significant and positive association between the mean yield and  S1, S2, S6, N1, N3, Gai, and Y as well as Pi_f, Pi_u, and Pi_a.It was demonstrated that there was a positive and significant correlation between the stability statistics N1 with S2, S6 with Pi_f, S6 with Pi_u, S6 with Pi_a, Pi_f with Pi_u, Pi_f with Pi_a, Pi_u with Pi_a, Pi_u with Gai, Pi_u with Y, Pi_a with Gai, and N4 with S1.

The multi-trait genotype-ideotype distance index (MGIDI) selection process and selection gains
The MGIDI function makes it easier to compute the index, allowing for identifying idiotypes by declaring which traits were desired to increase or decrease 48 .The lower contribution of a factor in the genotype indicated that the genotype is highly selective for the trait in a particular factor.The factors closer to the edge in Fig. 9 are less distant from the desired idiotype, and this can be used for genotype selection.
To identify the genotypes that conserve water, the multi-trait genotype-ideotype distance (MGIDI) index was calculated, accounting for all evaluated features.A significant genotypic influence was found for each of the five features that were measured: G1 (Binadhan-10), G2 (Binadhan-8), G3 (BRRI dhan-29), G4 (BRRI dhan-47), and G5 (Binadhan-17) (Fig. 9A,B).The analytical results showed that G3 (BRRI dhan-29) was selected for both The analysis involving two main factors from a set of six traits delineates their collective contribution to 58.4 percent of the total variance among traits, revealing specific trait clusters.These factors offer a simplified view of trait interconnections for AWD: Factor 1 includes traits such as ET, DM, and GY, while Factor 2 consists of PH, PL, and TGW for AWD.Communal and unique attributes account for 86.6% and 32.4%, respectively, of all the genetic diversity within the dataset, according to Table 5.This supports the effectiveness of factor analysis in creating an index that facilitates optimal trait selection 41 .The use of the MGIDI index resulted in favorable genetic enhancements across all examined traits, with an overall genetic gain of 13.06 percent.Notable traits like TGW and ET showed impressive selection gains of 4.78 and 5.55 percent, respectively (Table 6).This highlights the potency of the MGIDI index in promoting precise and advantageous trait selection, thereby improving crop development strategies.
The radar plot visualizations (Fig. 10A,B) provide a detailed analysis of the strengths and weaknesses of various rice genotypes across two moisture conditions.These figures are instrumental in evaluating the genotypes' adaptation by examining how each trait contributes to the Modified Genotypic Ideotype Deviation Index (MGIDI).Each factor's contribution towards the MGIDI is strategically ranked, with the most influential factors positioned closer to the center of the plot and the less influential ones farther from the center.Traits that are closer to the outer edge of the radar plot are considered more aligned with the ideotype.This implies that a smaller proportion of the genotype variation is explained by these factors, indicating a closer similarity to  www.nature.com/scientificreports/ the ideotypical traits.A view on strengths and weaknesses under AWD revealed that the performance of the selected genotypes, viz., G3 showed strengths related to factor FA1 that holds ET, DM, and GY, while Factor 2 consists of PH, PL, and TGW with positive SGs except TGW (Fig. 10A).This analysis provides a comprehensive understanding of how each genotype performs under varying moisture regimes and which traits contribute most effectively towards achieving an ideotype-like phenotype.This systematic approach aids in identifying key areas of strength and potential improvement for rice breeding programs under different environmental conditions.Similarly, under CF conditions (Fig. 10B), most of the selected genotypes that is G3 showed strengths related to all the factors (FA1 and FA2).According to FA1, genotype G3 performed very well with a maximum contribution towards the MGIDI.The examination of two primary factors from an initial group of six traits uncovered that they account for 68.8 percent of the variation among these traits, showcasing distinct clusters of attributes.This analysis narrows down the complexity of trait interrelations specifically for CF: Factor 1 (FA1) includes traits such as ET, DM, and GY, while Factor 2 comprises characteristics like PH, PL, and TGW in CF.The analysis displayed an average communality and uniqueness representing 95.3% and 5.6% of the genetic variance, as detailed in Table 6, affirming the utility of factor analysis in forming an index conducive to optimal trait selection 41 .Further, the implementation of the MGIDI index yielded favorable genetic improvements across the evaluated traits, achieving a total genetic gain of 16.52 percent.Particularly, traits like PL and ET marked significant selection gains of 8.52 and 3.42 percent, respectively (Table 6).This underscores the success of the MGIDI index in promoting specific and beneficial trait selection for advancing crop improvement strategies.

Discussions
For irrigated rice, the AWD irrigation system was created as a way to save water.The method has naturally been used mostly in places with limited water resources.In this study, water treatment had a dominant effect on grain yield, in addition to the effect of the cropping season, which was significant due to climatic changes.The International Rice Research Institute (IRRI) created the AWD technique, which was first presented in Bangladesh in 2004, to meet the demands of farmers who want to reduce their use of fuel, water, and energy in their irrigated rice production.Lowland rice producers may utilize the AWD, a water-saving technique, to use less water in their irrigated crops.AWD involves flooding the land with irrigation water a certain number of days after the ponded water evaporates.Water management at saturated soil conditions was shown to sustain similar moisture content in the soil, thus supporting leaf physiological and plant growth performances that resulted in maintaining a high grain yield of rice.In AWD, the interval between irrigations might range from one day to more than ten days during which the land is not watered."Safe AWD" refers to the water depth threshold of 15 cm (below the surface) before irrigation, as it prevents yield reduction.Water savings with Safe AWD range from 15.0 to 30.0 percent.The country's north-central, northwestern, and southwestern regions face even more pressing problems due to falling groundwater levels and growing water shortages.Thus, there is little doubt that the intense groundwater abstraction caused by dry season rice cultivation is responsible for the groundwater levels' decline, which is occurring at a rate of 0.1-0.5 m/yr 49 .The majority of morpho-physiological and agronomical features of rice cultivation are significantly impacted adversely by heat and/or drought, and this is reflected negatively in grain production, particularly if these stressors occur during crucial stages of the plant life cycle.Therefore, the creation of cultivars with high productivity together with high stability under fluctuating drought and/or heat levels is the primary focus of plant breeding specialists and the long-term goal of modern breeding programs 50,51 .Thus, in the present investigation, genotypes with high adaptability and stability are chosen using AMMI and GGE models.The environment and its interactions have a significant impact on yield and yield-related characteristics, therefore understanding the G × E interaction is crucial for identifying stable genotypes in breeding programs.
The yield component characteristics PH, ET, PL, TGW, and DM were very significant, which is consistent with the findings of Yadav et al. 52 .This suggests that each genotype is genetically diverse and reliable for character selection of these variables.The results of the ANOVA study showed that there were significant differences in the variances related to genotype, environment, and genotype × environment interaction for the characters PH, ET, PL, TGW, and DM.9][60][61][62][63][64][65] .The noteworthy interplay resulting from G × E serves as a gauge for the diverse reactions of rice genotypes to quantitative attributes in both scenarios.This indicates that to identify genotypes with either particular adaptation or broad adaptability, genotypes must be assessed for stability analysis 66 .The genetic variety of the genotypes, environmental factors, and the interactions between the genotypes and environments might all be attributed to the extremely significant variations in environments and genotypes found in this research.Numerous studies have shown comparable findings, and the variety of the multiple conditions examined with the various genotypes planted may account for the extremely substantial GEI for certain agronomic variables [67][68][69] .Numerous studies have shown that GEIs significantly affect rice yields 51,70,71 .
The AMMI model, which is the popular stability analysis method, measures the amount of variability in both genotype and environment using a combination of ANOVA and GEI 44,72 .The main effects (genotype and environment means) are shown by the X-coordinate in the AMMI 1 biplot (Fig. 2 A-D), while the interaction effects (PC1) are represented by the Y-coordinate.Significant differences exist across genotypes in the AMMI1 biplot concerning direction and amplitude along the X-axis (yield) and Y-axis (PC 1 scores).A larger score (absolute value) implies instability and is specific to particular situations when evaluating a biplot experiment.Main effects with a PC score around zero suggest low interaction effects and are stable 73 .A positive interaction occurs between a genotype and its environment when their signals on the PC axis match; a negative interaction occurs when their signs disagree 74 .Genotypes with a PC score of almost zero and a root average yield larger than the grand mean, according to the AMMI model, are often able to adapt to any environment.Nonetheless, it is believed that genotypes with high IPCA scores and mean performances are better suited to their surroundings 47 .AWD G3 showed particular flexibility for E7, E2, and E4 conditions with mean yields below the grand mean, as shown in Fig. 2A.The AMMI1 biplot explained the distribution of the tested genotypes based on the GY mean and PC1 values, allowing for the identification of the genotypes G3 and G2, which had the greatest GY mean and the lowest IPCA1 scores (Fig. 2B).The association between genotypes and environments and their environmental adaptability is shown by the PC score of genotypes 75 .The most stable genotypes are those with an IPCA score close to zero 76 .A genotype is more stable, according to Balakrishnan 77 , if its IPCA score is closer to zero.Conversely, a high IPCA score positive or negative indicates that the genotype has adapted to a certain environment.The AMMI2 biplot provided more explanation for the multiplicative effects of GEI throughout the first two IPCAs (Fig. 5C,D).As a result, the genotypes G3 and G5 were the best; they were found near the origin and had the lowest PCs values (Fig. 2C).Genotypes G5 and G1 showed practically zero scores for both PC1 and PC2, indicating that these genotypes were less interacting with the locations, according to the AMMI2 biplot (Fig. 5D) graphical analysis for PC1 and PC2.Our findings are in line with those of Gauch and Zobel 78 and Khan 69 , who discovered that the first two PCAs are enough for the AMMI model's projection, despite some research suggesting the use of the first four PCs in the multi-environment trial 79 .AMMI biplots were used in several research, such as those on barley 80 , rice 25 , sorghum 42 , and maize 74 , to identify the corresponding genotypes.
The average environment coordinate (AEC) technique was used to assess trait performance and genotype stability 46 .The average environment coordinate is the line that crosses the perpendicular line of the biplot origin.The direction in which the arrows point to a tiny circle represents the genotype mean yield 81 .In the meanwhile, genotype stability is shown by the perpendicular lines that cross the biplot origin 82 .Higher genotype stability is seen in genotypes that are nearer the AEC line 83 .The optimal genotype is stable and has a high mean yield per hectare.G3 and G5 genotypes were associated with a high mean yield per hectare (Fig. 5A).While G5 was unstable, as seen by longer lines from AEC, genotype G3 performed steadily with lines closer to AEC.While G1 was shown to be more closely connected to AEC lines, suggesting consistent performance and a better mean yield per ha production, Fig. 5B revealed that genotypes with a higher mean yield per ha were G3.The genotypes positioned in the inner circle were considerably more appropriate than those in the outside circle.The genotypes that are next to and closer to the inner circle are thought to be best; nevertheless, in a few cases, no genotype was discovered within the inner circle 69,84 .Figures 7A,B show the genotype ranking according to mean and stability routine.The AEC approaches were used in GGE biplot analysis to estimate genotype yield and stability 46 .The average environment coordinate (AEC) is the line that goes through the biplot origin.It is determined by averaging the PC1 (mean yield) and PC2 (stability) scores for every environment 45 .A greater mean yield is shown by being closer to the concentric circle.The best genotypes for selection are those that have a high mean yield and good stability.They have the smaller vector from the AEC and are situated near the origin in the biplot.The genotypes with yield performance larger than the mean yield are on the right side of the line without an arrow, whereas the genotypes with yield performance less than the mean yield are on the left side of this line.The most stable and highly productive genotypes in this investigation were G3 in both treatments (Fig. 6A,B).Each environment's rating is shown by how close the environmental point is to the ideal point (small arrow) in Fig. 7 (environment rank).Ruswandi 85 states that the perfect environment is the one that is closest to the optimum point.The location that is closest to the ideal point, E2, was the environment that was outside the circle and had a sufficiently wide angle with the ideal point in our study of the AWD.This was followed by E7 and E1 (As seen in Fig. 7A).This suggests that although other settings may be utilized to select adaptable genotypes, the E2 was the most suitable environment for selecting stable and great-producing genotypes.In Fig. 7B, the E9 sites are the most wanted because they are the closest to the ideal environment, while the E7 locations were thought to be the most undesired.
According to Azam et al. 86 and Hossain et al. 87 the what-won-where perspective of the GGE bi-plot is the best model for multi-environment trial data.It enables you to classify the environments and identify the genotype that performs best in each.This biplot piqued the interest of many academics since it presents key ideas including crossover GE, massive environment difference, and specialized adaptation in a visual manner.To include all additional genotypes, a polygon is initially created on the genotypes that are furthest from the biplot origin.Beginning at the biplot origin, draw the perpendicular lines to each side of the polygon 88 .In the first mega-environment with AWD, Vertex genotype G3 emerged as the dominant genotype.However, in the second mega-environment, genotype G5 emerged victorious (Fig. 4A).In the mega-environment of CF therapy, genotype G3 fared better as well (Fig. 4B)."Different cultivars should be selected and deployed for each different environment," the conclusion states.The rice workers, namely Akter 82 , Mohan 89 , and Siddi 90 , reported similar findings.
The MGIDI is a relatively new technique that is being used in an increasing number of crops to assist in discovering a better genotype for many characteristics simultaneously 86,87,91,92 .The Euclidean genotype-ideotype distance is computed using a factor analysis scores 93 .The MTSI is based on the projection of factor analysis results onto the genotype-ideotype distance.Plant breeders may find the MGIDI beneficial in choosing genotypes based on numerous characteristics since it offers a robust and understandable selection strategy.Based on the characteristics under evaluation, the genotypes with lower MGIDI values have stronger stability.Six characteristics were taken into account for MGIDI analysis in our investigation: G3 genotypes were thought to be the most desirable or stable.The MGIDI index is a unique and easy-to-interpret selection procedure that has many practical applications to obtain long-term genetic gain in primary traits (such as grain yield) without jeopardizing gains of secondary traits (such as plant height).The proportion explained by each factor towards the MGIDI index, i.e., "strengths and weaknesses view" (Fig. 10A,B), is an important graphical tool to identify the strengths and weaknesses of test hybrids in terms of "trait (group of traits) need to be improved" and is an added advantage over existing indices.A similar study was reported by Olivoto and Nardino 41 , who evaluated a set of 13 strawberry cultivars, wherein the strengths were described by employing MGIDI.

Description of the site and soil
In Bangladesh, the Boro season lowers the water table and calls for greater irrigation.The shortage of water in Bangladesh's western Borendra area makes rice production more difficult.The ideal sites for this experiment were identified using surface reflectance (SR) data obtained from the Moderate Resolution Imaging Spectra Radiometer (MODIS) satellite equipped with a TERRA sensor and a product called MOD09.A sinusoidal grid with a horizontal resolution of 250 m is used to protect all of the data.The data was collected using the web-based warehouse inventory search tool WIST between 2010 and 2019.Month-by-month maps with the relevant photos were created when the data analysis was finished (Table 7 and Fig. 11).
The High Barind Tract (HBT), which covers an area of around 1600 km 2 (159,964 hectares), is located in Bangladesh's North-West Rajshahi Division at coordinates of 24.40-24.80°N and 88.30-88.80°E. Its terraced terrain, poor fertility soil, sparse vegetation, lack of main river routes, and relatively low rainfall combined with a protracted dry season from October to May set it apart from other regions of the nation 94 .Because of its high summertime air temperatures, which can reach maximums of 35-40 °C in April and May, its limited groundwater reserves and recharge, the poor water-holding capacity of surface soil during the post-rainy season, and its relatively low and irregular precipitation (average annual rainfall of 1440 mm, which can vary spatially and seasonally between 790 and 2200 mm), it is known as the most drought-prone region in the country (Fig. 12).Unlike Bangladesh's floodplain regions, which experience yearly flooding, the HBT is comparatively elevated (9-47 m above sea level 94 .
The High Barind Tract's soils reacted with a mild acidity.According to the USDA classification, they were of medium-to heavy texture, with textural grades of silt loam for six soils, silty clay loam for three soils, and silty clay for the remaining soil.The majority of soils have significant clay contents, ranging from 16.0 to 43.9%, with topography being the primary cause of the difference.Quartz predominated in the 2-20 gm silt fraction, with minor amounts of mica, plagioclase, K-feldspar, and chlorite 95 .Generally speaking, the texture of the soil was clay to silty clay underneath and silty clay to clay on top.In every profile, there was a gradual rise in the amount of clay from the surface down 96 .

Plant components and the layout of field experiments
The present investigation used AWD and CF systems, whereby the five genotypes including the genetic materials better than the others (Table 8).Table 8 displays the genotype name, code, collection source, and varietal description.The field experiment was performed from January to May 2020, throughout the Boro season.30 days after the seedlings were grown into the seedbed, they were moved per hill, with a 20 cm × 20 cm gap between rows and plants 97 .Plants were spaced 20 cm apart from one another in this experiment, and rows were spaced 20 cm apart.The nine settings in the multi-environment trials were performed under alternating cultivations of early and standard planting seasons.Rice cultivars and irrigation schedules were the two variables combined, and they were set up in a split-plot design with three replications.The main plot was the irrigation regime, while the subplot was the rice genotype.Rice cultivars were assigned at random to three replications at a subplot.A 6-m buffer was placed between the two irrigation regime plots to prevent sub flow from the continuous flooding plot from contaminating the plot.The experiment used a split-plot design with three (3) replications, maintaining a unit plot size of 8 × 5 m.To keep the soil fertility in the experimental units at the optimal level for rice cultivation in Bangladesh, additional chemical fertilizers were sprayed at rates of 220-120-90-60 kg/ha, respectively, including urea, triple super phosphate, muriate of potash, and gypsum (CaSO 4 .2H 2 O) 98 .The full doses of TSP and muriate of potash, along with half of the dose of urea, were applied as the basal fertilizer application, and the remaining 50% of the nitrogenous fertilizer was divided into two halves.The first one was used during tillering, and the second one was used during booting.All crucial intercultural tasks, such as gap filling, routine weeding, applying insecticides, adding fertilizer to the top, and so on, were completed throughout the growing season to support the plants' healthy growth and development.Experimental research and field studies on rice including the collection of rice seeds comply with our institutional, national, and international guidelines and legislation.
Two AWD irrigation regimes were used as a water-saving strategy in place of constant flooding.Throughout the whole cropping season, under the continuous flooding regime, a constant water layer of 5 to 10 cm was maintained, depending on the crop growth stage (Fig. 13a).Irrigation was the first water-saving measure considered when the soil matric potential (AWD30) reached 30 kPa (Fig. 13b).The average daily soil matric potential (SMP) was measured and hourly monitored using Watermark Granular Matrix sensors (WGMs; model: EC06-783, Irrometer, Co., Riverside, CA, USA).WGMs provide an indirect way to quantify SMP by directly measuring soil water tension.Two WGMs were placed into each plot at a depth of 15 cm to lessen the variation in soil water content caused by field leveling.As soon as the standing water level hit 20 mm, irrigation was turned on while using the continuous flooding irrigation technique.Thirty millimeters of irrigation water were applied to the rice during its vegetative growth stage, and eighty millimeters were applied during its reproductive phase.During the rice vegetative stage, 30 mm of irrigation water was applied at each irrigation event under the AWD treatments (AWD30).During the rice reproductive phase, irrigation was similarly managed under AWD treatments to continuous flooding.This was done to avoid or reduce sterility.Before imposing irrigation regimes on the treatments, 100 mm of irrigation water was administered for paddling and transplant survival 99 .

Data gathering
In each plot, 10 randomly selected plants were taken for sampling all of the phenotypic data.Plants were harvested when 85% of the grains of the panicles turned a golden yellow color.Three biological replicates (10 plants per replicate for each trait) were used to collect data on growth and yield-related traits.The studied days to maturity (DM) trait was measured from the sowing to harvesting day.At harvest, plant height (PH) was determined from the soil surface to the highest part of the plant; Effective tiller (ET) was counted as the number of panicles bearing tillers per hill; Panicle length (PL) was determined from the bottom of the peduncle to the topmost of the Table 7. List of the locations selected for the experiment location detail during October 2019 to April 2020.(Source: Bangladesh Meteorological Department (BMD), Rajshahi region).panicle; Thousand seeds weight was (TSW) recorded by counting one thousand seeds manually from each pot and weighed using an electronic balance.The grain yield was obtained using a crop cutting test (CCT) of 1 m by 1 m areas of each experimental plot.The harvested grains were dried, winnowed, and weighed.The weight was then converted to per unit area crop yield based on 14% grain moisture content and presented as grain yield (t/ ha).Tons per hectare were calculated from the seed yield of each plot.

Statistical analysis
The average data of various traits (Supplimentary file 1 and 2) were analyzed statistically 100 and biometrically 101,102 .R version 4.2.1wasused for all statistical analysis and data visualization 103 .The analysis of variance, AMMI analysis of variance 104 , genotype plus genotype × environment (GGE) biplot 46 , and stability statistical analysis were all performed using the "metan" R Package 93 .The AMMI model uses two techniques: multiplicative interaction effect and additive main effects of environments and genotypes, as well as ANOVA and singular value decomposition 105 .
A genotype main effect plus genotype-by-environment interaction (GGE) biplot was used to display the winning genotypes and mega-environments (MEs) for yield-related features and grain production 46,106 .Using a "which-won-where" biplot, investigations on MEs and the optimal genotype in each habitat were carried out."Discriminativeness vs. representativeness" biplots were used to assess the ideal conditions for each characteristic.Ideal genotypes exhibiting high mean performance and stability across all parameters were shown using

Figure 8 .
Figure 8. Heatmap-based Speraman's correlation coefficient among various stability parameters with yield data of 5 rice genotypes evaluated at 9 test environments.Parametric measures such as: ASV-AMMI stability value, RPGV: Relative performance of genotypes values; HMRPGV: Harmonic mean of RPGV, Ecoval: Wricke's ecovalence.Non-parametric measures such as: N1, N2, N3, N4: Thennarasu"s statistics, Sij: Deviations from the joint-regression analysis, S1: mean of the absolute rank differences of a genotype over the n environments, S2: variance among the ranks over the k environments, S3: sum of the absolute deviations, S6: relative sum of squares of rank for each genotype, Gai: Geometric adaptability index, Pi_a, Pi_f, Pi_u: Superiority indexes for all, favorable and unfavorable environments, respectively.

Figure 10 .
Figure 10.The strengths and weaknesses view of the selected rice genotype is shown as the proportion of each factor on the computed MGIDI values over two moisture regimes, namely, (A) AWD and (B) CF.

Figure 11 .
Figure 11.Water status of soil from October to April in Bangladesh.The data were collected from 2010 to 2019 by using the web-based warehouse inventory search tool (WIST).

Figure 12 .
Figure 12.Temperature, humidity, and rainfall variation during grain filling period in both AWD and CF methods at 9 environments in the rice crop season of 2019-2020.

Figure 13 .
Figure 13.Comparison of canopy of crop and soil wetness status under two different irrigation regimes i.e., continuous flooding and alternate wet and dry conditions during the rice tillering stage at Nachole during Boro season.(A) Continuous flooding, (B) Alternate wet and dry.

Table 5 .
Eigenvalues, explained variance, factorial loadings after varimax rotation, and communalities obtained in the factor analysis within the soil moisture regime based on MGIDI values for AWD.

Table 6 .
Eigenvalues, explained variance, factorial loadings after varimax rotation, and communalities obtained in the factor analysis within the soil moisture regime based on MGIDI values for CF.

Table 8 .
List of the Genotypes selected for this experiment.